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Abstract 

Experimental data on the pp — > pmr~^ reaction measured in an exclusive two- 
arm experiment at 800 MeV show a narrow peak arising from the strong proton- 
neutron final-state interaction. It was claimed, within the framework of a certain 
model, that this peak contained up to a 25% spin-singlet final state contribution. 
By comparing the data with those of pp dvr"'" in a largely model- independent 
way, it is here demonstrated that at all the angles measured the whole of the peak 
could be explained as being due to spin-triplet final states, with the spin-singlet 
being at most a few percent. Good qualitative agreement with the measured proton 
analysing power is also found within this approach. 
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Inclusive data on the pp n^X reaction in the 400 MeV to 1 GeV range 
generally show two structures. There is firstly a broad peak corresponding to 
the quasi-free production of the A-isobar. In addition, there is an enhancement 
near the edge of phase space arising from the strong neutron-proton final- 
state interaction (fsi) in either the spin-triplet or singlet 5'-wave which comes 
when the np excitation energy, En,,, is only a few MeV. The details of the 
enhancement region are hard to investigate in a single-arm experiment, where 
only the tt"*" is measured, because of contamination from the much larger two- 
body — > dTT"*" reaction [1]. The simplest way to overcome this background 
is by measuring in coincidence the final proton and pion in the exclusive 
pp — ^ pn 7r+ reaction. The most complete examples of such an experiment were 
carried out at LAMPF, where both the five-fold differential cross section [2] 
and the proton analysing power [3] were studied at a proton beam energy of 
Tp = 800 MeV. 

LAMPF data taken at one pair of proton/pion angles are shown in Fig. 1 
as a function of the detected proton momentum [2]. The large peak on the 
right is associated with protons and pions formed in the pp — > A+^n/A+p 
reactions and its magnitude and width can be explained in different versions 
of one-meson-exchange models [2,4]. At the maximum of the smaller peak, 
Enp is about 1 MeV and the cross section is strongly influenced by the np fsi. 
The form of this peak is normally parameterised by Watson final-state inter- 
action factors [5,6] which take into account the nearby poles in the scattering 
amplitudes due to the deuteron bound state in the spin-triplet case and the 
anti-bound state for the spin-singlet. Since the latter is closer to E^p — 0, the 
singlet fsi peak is expected to be narrower than that of the triplet. Though this 
difference in shape can, in principle, be used to extract the relative amount of 
the np spin-triplet and singlet in the final-state peak, the limited number of 
points in the peak and the modest resolution in E^p makes this impractical 
for the LAMPF data. 

In analysing their experiment, the authors of Ref. [2] made the ad hoc assump- 
tion that the one-meson-exchange model for quasi-free A production was valid 
for Enp above 10 MeV and that, below this excitation energy, the prediction 
could be smoothly joined onto Watson fsi factors. The relative spin-singlet 
strength depends, of course, upon the kinematics and their analysis suggested 
that the fraction was about 15-25% of the total for the different angle pairs 
measured. This is close to a statistical mixture of 25% and is in complete con- 
trast to measurements at lower energies, where the singlet fraction is typically 
10% [7]. Alternative model-dependent analyses of the LAMPF data do lead to 
smaller singlet production, but they depend upon other assumptions made [8]. 
In view of these differing conclusions, it is worthwhile to seek a different way 
of deducing the spin-singlet contribution from these data. 

When the square of the low energy np triplet 5'-state scattering wave func- 
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tion at energy Enp = k'^/niN, where mjv is twice the n-p reduced mass, is 
analytically continued to negative energy, it manifests the deuteron pole at 
— — q;|, where at = 0.232 fm~^. It has recently been shown [9] that the 
relative normalisation of the scattering, ^'^(r), and bound-state, \E'a((r), wave 
function depends purely upon the deuteron binding energy and is independent 
of the np potential [9]. Using real boundary conditions, it follows that at short 
distances 

[^k{r)f^ [^a.(r)f ■ (1) 

Q;t(/c^ + af) 

Apart from questions associated with the D-state contributions which are 
small at low /c^, this relation becomes exact as /c^ — > — af. In the scattering 

region where Enp > 0, it allow one to estimate the spin-triplet amplitude for 
pp pmi^ in terms of that for pp dn^. Thus, for any value of the initial 
(aj) and final spin-triplet projection A, the matrix elements are related by 

Ka^iPP ^ {nphn-") ^ f{k') M^^^^ipp ^ dn+) , (2) 

where 

ne) = • (3) 

We are here using a normalisation where the unpolarised pp — > dn'^ cross 
section is given by 

^. = ^'f\T.\MUw^d.-)\\ (4) 

USi U^TT App (/p ^ ^^^^^ 

Spp is the square of the centre-of-mass (cm) energy, and g* and g*^ are the cm 
momenta in the initial and final states respectively. At low Enp the approxi- 
mation of eq. (2) reproduces very well the results of single-arm experiments 
where only the pion is detected [7] . 

We now extend this approach to treat two-arm experiments. Starting from 
eq. (2), the triplet contribution to the laboratory five- fold differential cross 
section for the detection of a pion at an angle and a proton at an angle 6p 
becomes 

dp, dO. - 167r3 '^ql, * ^ ^ dQ* ^ ^ ' 
where the phase-space factor is 

$ = ^ . (Q\ 

PO rUp Ep \pI En-Pn- PnEn \ 
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Here po is the beam momentum, Ei and pi {i = p, n, n) are the laboratory 
energy and momentum of the i-th particle in the final state. 

The pion production angles in the laboratory (9^^) and cm (9*) systems are 
related by 

Eo - Po p^ cos 9^ =Soe^- q* ql cos 9^ , (7) 

where £7r(£o) ^-nd g* are the energy of the pion (incident proton) and 3- 
momentum of the pion in the overall cm system and Eq is the total laboratory 
energy of the incident proton. Due to the difference between the effective mass 
of the final pn state and that of the deuteron, ^ q*^. However, this effect 
is quite small at low Enp- 

The results of eq. (5), which should be valid at low relative energies Enp, do not 
depend upon the details of the pion production dynamics and automatically 
include the final state interaction in the triplet pn system. 

The input pp ^ dir^ cross sections are taken from the SAID SP96 parameteri- 
sation [10] and this procedure should involve errors that arc smaller than those 
of an individual experiment. The predictions of our approach for the triplet 
VP ~^ cross section are shown at one angle pair in Fig. 1 together with 

an evaluation of quasi- free A-production [4]. The values of the np excitation 
energies are indicated and from this it is seen that when Enp is below about 
5 MeV the magnitude and shape of the fsi peak are both well reproduced. For 
comparison the fsi parameterisation of Goldberger and Watson [6] , multiplied 
by phase space and normalised to the peak value, is also illustrated. 

Results at different angles in the fsi region are shown in Fig. 2. Larger devi- 
ations from the predictions are to be found when the minimum value of the 
excitation energy, -E™", is increased. This is the case at {9p,9.„) = (25°, 40°) 
and (30°, 28°) where are respectively 13 and 28 MeV. Even in this last 
case, the changes induced by the modified kinematics are quite small, as can 
be judged from the dot-dash curve shown in the figure. A potentially more se- 
rious effect arises from averaging the estimates over the experimental angular 
acceptance. As shown in Fig. 2, this tends to reduce a httle the predictions 
when -B^*" < 1 MeV but can increase them otherwise. After smearing, there 
seems to be some undcrprediction in Fig. 2b, though it should be stressed 
that these data were taken from the polarisation experiment [3] where the 
consistency with the earlier unpolarised cross section run was found only on 
the 15% level. 

Since the predictions of the Faldt-Wilkin extrapolation theorem [9] reproduce 
most of the magnitude and angular dependence of the cross section leading to 
the fsi peaks, this is strong evidence that the vast bulk of the LAMPF data 
corresponds to np triplet final states. Extra confirmation of this interpretation 
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is found from the proton analysing power which was measured at two pairs 
of angles in the fsi region [3]. In the fsi peak we expect Ay to be essentially 
constant at the value corresponding to that of pp — > d tt"*" at the appropriate 
pion angle. These values are shown with the experimental data in Fig. 3. 
Though at (14.5°, 21°) the data seem to fall a little above the prediction, 
deviations of this size are not unknown in polarisation measurements. 

Although we have reproduced well the LAMPF data at low Enp with just 
triplet terms, it may be helpful to try to estimate an upper bound on the 
possible spin-singlet contribution from the areas under the peaks. The overall 
systematic error due to beam normalisation and detector efficiency in the mea- 
surement was estimated to be about 7% [2], to which must be compounded 
some error coming from the pp ^ d7r~^ input [10]. The error arising from us- 
ing the extrapolation theorem in the scattering domain is likely to be rather 
smaller than this, with variations in the wave functions at short distances 
being on the 1-2% level for Enp < 10 MeV [11]. Further work is needed to 
include the D-state effects more consistently, though it has been shown that 
the extrapolation theorem is valid in this coupled-channel case [12]. Under 
the present conditions, we would estimate the accuracy of eq. (5) to be better 
than 5% for E^p < 3 MeV [13]. We therefore conclude that the singlet contri- 
bution to the unpolarised LAMPF data at low Enp in Fig. 1 is below about 
10%. This upper bound is about a factor of two less than the average quoted 
in the experimental paper [2], though it must be stressed that this involved 
considerable model dependence, including the choice of a 10 MeV matching 
point. 

There is a 15% normalisation uncertainty between the first measurements of 
the unpolarised cross section at LAMPF [2] and their later experiment, where 
the primary purpose was the determination of the analysing power [3]. If, 
nevertheless, we take seriously the disagreement with the latter data shown in 
fig. 2b, it is possible that this is due to a spin-singlet contribution. Assuming 
that the singlet cross section to be of the form 



where = —0.04 fm~^ then, after averaging over acceptance, a value of 
C = 0.03 restores the agreement with the data. This would correspond to a 
singlet fraction integrated over the range < Enp < 3 MeV of about 10%. 
However, such a fraction would make the agreement worse in fig. 2c, where 
the data were taken at a rather similar cm angle but in the other hemisphere. 

To improve the sensitivity to the singlet/triplet ratio using the extrapolation 
theorem, the ratio of the pp — > pmr^ and pp ^ dir^ cross sections has to 
be established better by measuring both reactions in the same experiment. 
Data with better resolution on Enp would also allow one to investigate the 
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singlet/triplet ratio from the shape of the fsi peak. A new measurement of the 
pp — > pn7r~^ reaction at the ANKE spectrometer of the proton synchrotron 
COSY-Jiihch, where both protons and pions were detected near the forward 
direction at Tp — 492 MeV, involved both these improvements and can there- 
fore put more stringent hmits on the spin ratio [14]. 

In conclusion, we have generalised the use of the extrapolation theorem linking 
np scattering and bound-state wave functions to describe five-fold differential 
cross sections. In so doing, we have shown that the LAMPF pp pn 7r+ cross 
section data at 800 MeV are consistent with there being no final-state spin- 
singlet contribution. This smallness should not come as a complete surprise. 
If one assumes that pion production passes through an intermediate A iso- 
bar [2,4,8], then the 5'-wave A intermediate state cannot lead to singlet np 
final states. These must be produced through higher partial waves or from 
non-resonant pion production. 

This work was initiated within the framework of the ANKE collaboration at 
COSY and the authors would like to thank all their colleagues, especially 
V. Abaev, V.I. Komarov, V. Koptev, and V. Kurbatov, for many fruitful dis- 
cussions. Financial support from the Forschungszentrum Jiilich and BMBF 
(WTZ grant KAZ 99/001) and the generous hospitality of Institut fiir Kern- 
physik, where the work was carried out, are gratefully acknowledged. 
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Fig. 1. Differential cross section for the pp pmr~^ reaction at 800 MeV at fixed 
proton and pion laboratory angles of {6p, 9.,^) = (15°, 20.8°) as a function of the mea- 
sured proton momentum [2] . The pp AN peak on the right can be described in 
a one-meson-exchange model (dashed curve) [4] . The peak on the left is a reflection 
of the strong np fsi, which can be calculated (solid curve) from the pp (ivr"*" cross 
section, as discussed in the text. A scale is given showing the np excitation energy 
at the geometry corresponding to the centres of the counters. For comparison, the 
triplet fsi factor in the Goldberger and Watson form [6], multiplied by the phase 
space factor of eq. (6) and normalised to the peak, is shown as the dotted curve. 
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Fig. 2. Differential cross section for the pp pnir^ reaction at 800 MeV 

in the neighbourhood of the np fsi region for different angular positions, 
a) {ep,e^) = (15°, 20.8°), b) (14.5°, 42°), c) (14.5°, 21°), d) (25°, 40°), e) (20°, 22°), 
f) (30°, 28°). The np excitation energies at the counter centres are indicated; the 
number with the decimal corresponds to the minimum value of Enp. The experi- 
mental data are taken from Refs. [2] (circles) and [3] (triangles) and it should be 
noted that these may differ in normalisation by 15%. The data are compared to the 
spin-triplet final-state predictions of eq. (5) without acceptance corrections (solid 
curve) and with (broken curve). An attempt to resolve the apparent discrepancy 
in b), through the introduction of a singlet contribution by eq. (8) with C, = 0.03 
makes the situation worse in c) (dotted curves). The dot-dash curve in f) was ob- 
tained using the deuteron mass rather than the invariant mass of the np system 
when evaluating the pion production angle ^* in eq. (7). 
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Fig. 3. Proton analysing power for the pp pmr~^ reaction at 800 MeV in the 
neighbourhood of the np fsi region. The experimental data of Refs. [3] are compared 
to the SAID predictions of Ay (horizontal lines) for the pp ^ dir^ reaction [10] . It 
should be noted that SAID uses the opposite convention for the sign of Ay. 
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